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TRAVERSING THE TRAILING VORTICES OF HEAVY 
TRANSPORT AIRPLANES 
By William A. McGowan 

SUMMARY 


Results are presented of normal-load-factor calculations made for 
a light normal-category airplane and a light transport-category airplane 
traversing the trailing vortices generated by each of three heavy trans- 
port airplanes. With each light airplane, the normal load factors were 
determined for several penetration paths lying in a plane perpendicular 
to the trailing vortices and for three center-of-gravity locations and 
velocities. Also determined for the light normal- category airplane were 
the elevator deflection required to maintain 1 g flight and the vertical 
displacements of the airplane from the prescribed penetration paths while 
traversing the vortices . 


The methods used (formulated from available theories) and an example 
illustrating use of derived charts for computing normal load factors are 
given in the appendixes . 


The results indicate that light airplanes traversing the wakes of 
currently operational heavy transport airplanes can experience loading 
conditions that exceed the design limit and, in some cases, the design 
ultimate load factors. For light airplanes traversing the wake generated 
by a proposed supersonic transport airplane the design ultimate load 
factors can be greatly exceeded. It was also shown that load factors 
imposed by the vortex system could be alleviated by elevator deflection; 
however, because of inherent pilot and control-system lag, the load 
factors would almost certainly be aggravated rather than alleviated. 


INTRODUCTION 


The adverse effects on an airplane traversing the wake of another 
airplane are well known (see, for example, refs. 1 to 4) and a number of 
papers have been written which bear on various phases of the problem 
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(for example, refs. 5 to 8). As a result of these works, it is known 
that the vortex trails can persist for a minute or more under calm 
atmospheric conditions and are rapidly broken up by turbulence and that 
there is good agreement between the measured and the calculated intensity 
of trailing vortices. Also, these investigations reveal certain trends 
from which it can be reasoned qualitatively that (l) since the vortex 
trails are invisible no appropriate evasive action can be taken by the 
pilot except to be alert on calm days in high- intensity traffic areas 
whether or not another airplane is in sight, (2) reducing speed or flying 
either above or below the path of an airplane should reduce the magni- 
tude of the load factors, (}) either entering perpendicular to or 
quartering the vortex trails will result in approximately the same magni- 
tude of load factors except as such maneuvers might cause the period of 
the vortex impulses to coincide more nearly with the natural oscillations 
of the airplane, and ( 4 ) the effect of elevator motion on the load factors 
could either be canceling or aggravating, depending upon the phasing of 
the elevator motion with the impulses from the vortices. 

Although these factors and trends are known, very little informa- 
tion is available on either the actual quantities involved or a reasonably 
simple method for computing the load factors in a given case. Therefore, 
it is the purpose of this report to present the results of a series of 
systematic normal- load- factor calculations and to present the methods 
used. The calculations were made for a light normal-category airplane 
and a light transport-category airplane traversing, under various condi- 
tions, the wakes generated by each of three heavy transport airplanes. 

A list of the symbols used throughout the paper is given in 
appendix A. 


SCOPE OF CALCULATIONS 


The variables involved in this normal- load- factor investigation are 
contingent upon (l) the characteristics and flight conditions of the 
generating airplanes and ( 2 ) the velocity, static stability, penetration 
path, configuration, and elevator motion of the penetrating airplanes. 

The three heavy transport airplanes generating the wakes are referred 
to as the generating airplanes. The two light airplanes penetrating the 
wakes of the heavy airplanes are referred to as the penetrating airplanes . 
Three-view sketches of the penetrating airplanes are given in figure 1. 
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Generating Airplanes 

The characteristics of the generating airplanes required to establish, 
by theoretical means, the strength and behavior of the wing-tip trailing 
vortices are given in table I. Two of the three transport airplanes 
selected, a military transport Ga and a swept-wing civil transport Gg, 
are presently operational and the third airplane Gc represents a proposed 
supersonic transport of the future (ref. 9). Three generating airplanes 
were selected so that the relative effects of airplane geometry and weight 
on the normal load factors experienced by a penetrating airplane could be 
determined. The flight conditions are those that might exist shortly 
before landing or after take-off, as the probability of one airplane's 
traversing the wake of another would be greatest near airfields where 
the relative intensity of airplane traffic is high. 

A schematic sketch of the wing-tip rolled-up trailing vortices shed 
by the generating airplanes is shown in figure 2. The generating air- 
planes were considered to be flying straight and level for the present 
analysis . 


Penetrating Airplanes 

Airplane Pa is a typical light normal-category airplane of about 
2,000 pounds with externally braced high wing and single engine. Air- 
plane Pb is a light transport- category airplane in current use. (See 
fig. 1.) The geometric and aerodynamic characteristics pertinent to the 
analysis are listed in table II. The aerodynamic characteristics were 
obtained from either available methods or suggested values presented in 
the literature. 

Airplane Pa was considered to meet standards for normal-category 
aircraft as defined in reference 10. Airplanes in this category are 
intended for nonacrobatic operation. The positive and negative design 
limit load factors are 3*8 and -1.5; respectively. The limit load 
factors represent the maximum loads anticipated in service. The air- 
plane structure is, however, designed to support ultimate loads of magni- 
tudes li times the limit loads. The positive and negative ultimate load 
factors are, therefore, 5-7 and- -2.25, respectively. 

Airplane Pg, in the transport category, was considered to comply 
with the requirements of reference 11. The positive and negative design 
limit load factors are 2.5 and -1.0, respectively, and the corresponding 
ultimate load factors are 3*75 and -1.5* 
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The penetrating airplanes, prior to entering the trailing vortices 
of the generating airplanes, are trimmed for flight along the prescribed 
penetration paths. All the penetration paths are in the plane shown in 
figure 2 . The plane of the paths is perpendicular to the trailing vor- 
tices of the generating airplane. The line of vortex centers, noted in 
figure 2, is the reference line used in designating the level and inclina- 
tion of the penetration paths. An illustrative sketch of the various 
paths through the vortices is shown in figure 3 . 

Although the altitude along the penetration paths varies, the mass 
density of air is assumed to be constant at a value of 0.002309 slug per 
cubic foot, occurring at the pressure altitude of 1,000 feet. For con- 
sistency of calculations, the penetrating airplanes enter the vortices 
from the left, as indicated in figure 3* The flight conditions for the 
penetrating airplanes are enumerated in table III as cases 1 to 16. 

The effects of penetrating-airplane velocity, center-of -gravity 
location (static stability), penetration path, configuration, and vortex 
intensity on the normal load factor were determined from the load-factor 
time histories calculated for cases 1 to 14 by varying the aforementioned 
variables in a systematic manner. The elevator of the penetrating air- 
plane was held fixed during evaluation of these effects. 

The additional influence of an estimated elevator motion on the 
load factor of a penetrating airplane was determined as case 15* The 
elevator motion represented an attempt by the pilot to alleviate the 
normal load factors imposed by the vortices. Corrective pitch action 
(down elevator) was initiated as the penetrating airplane experienced 
the positive load factors upon first entering the vortex system. The 
estimated control motion was made consistent with what were believed to 
be reasonable control rates and pilot and control-system lag. 

In addition, calculations were made to evaluate the incremental 
elevator motion required to maintain 1 g flight while traversing the 
vortices (case 16). Also, vertical displacements of the airplane from 
a prescribed penetration path through the vortices were determined for 
cases 1, 2, and 3- 


METHODS AND RESULTS 


As defined in this paper, the normal load factor on an airplane 
traversing trailing vortices is composed of the increments related to 
steady-state flight, vortex flow, and airplane dynamic response in pitch 
due to the vortex flow and/or incremental elevator deflection. The air- 
plane is considered to traverse the vortices in such a way that the com- 
ponent loads induced are symmetrical about its center line and its 
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response is in the vertical plane only. The methods for calculating the 
vortex strength and normal load factors make use of various established 
theories. Details of the methods used are presented in appendix B. 


Trailing Vortices 

In order to evaluate the normal-load- factor increments related to 
the vortex flow and the airplane pitch response (motivated by the vortex 
flow), it is first necessary to determine the intensity and behavior of 
the vortices . 

In accordance with theoretical expressions given in reference 12, 
and in conjunction with the geometric characteristics and flight condi- 
tions of the generating airplanes, the vortex parameters were evaluated 
for a distance e behind the wing at which the wing-tip trailing vor- 
tices became essentially rolled up. These parameters are: (l) the 

vortex-core radius r c , (2) the tangential velocity w v of the circula- 
tory flow in the plane of the penetration paths (example shown in fig. 4), 
(3) the distance d between the core centers, (4) the distance x 
from the nearest vortex core where the airplane is considered to first 
penetrate the vortices, and ( 5 ) the downward velocity of the vortex 
cores w& with reference to the generating airplane's flight path. 

These vortex parameters define the intensity, velocity distribution, 
location, and movement of the trailing vortices for the generating air- 
planes and are given in table IV. 

The theoretical expressions indicate that the tangential velocity 
approaches infinity near the vortex center. However, in this paper the 
vortex core is assumed to rotate as a solid body; hence, the velocity 
within the core decreases linearly from a maximum at the periphery to 
zero at the center. 


Normal Load Factor 

The normal load factor on the penetrating airplanes traversing the 
vortex fields can be written 


n = n x + n z + n r (l) 

where nq represents the trim steady-state or datum load factor, nj 

represents the load factor related to the vortex-system velocity, and 
n r represents the load factor due to airplane motion or dynamic response 

about the penetration path. Equation (l) can be written in coefficient 
form as 


1 . 
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n * 1 + (sr f + °* ) {2) 

The vortex- induced angle of attack on the penetrating airplane a i was 
found by determining the ratio (angle) of the resultant vortex-system 
velocity component (perpendicular to the penetration path) to the forward 
velocity of the penetrating airplane and then adjusting this ratio (angle) 
for lag-in-lift effect. An example of the resultant vortex-system velocity 
component for case 2 is shown in figure 5 • The induced angle of attack 
Op due to airplane response about the penetration path was evaluated 
from the second-order differential equation for the longitudinal motion 
of an assumed rigid penetrating airplane with impressed moment caused by 
the vortex flow and/or elevator motion. The values of the other param- 
eters in equation (2) are listed in table II. Equation (2) was solved 
for n at times selected along the penetration paths to represent ade- 
quately the normal-load- factor time histories. Typical normal load 
factors nq, nj, and n r and the sum of these factors calculated for 
an elevator-fixed penetration case are shown in figure 6. 

Plotted time-history results are shown in figure 7 for the calculated 
normal load factors on airplane Pa traversing with elevator fixed the vor- 
tices of airplane Ga- The time histories are grouped to illustrate the 
effects of velocity (fig. 7(a) ), center-of -gravity location (fig. 7(t)), 
height of penetration path with respect to the line of vortex centers 
(fig. 7(c)), and angle of penetration with respect to the line of vortex 
centers (fig. 7(d)). In figures 8 and 9 the normal-load-factor time 
histories are presented for airplanes Pa and P^ ; respectively, when 
penetrating with elevator fixed the vortices of airplanes Ga, Gb, and Gc • 
The effect of an estimated elevator motion on penetrating airplane Pa is 
shown in figure 10. Also shown in the figure is the time history of the 
estimated elevator motion. In addition, the time history of the initial 
incremental elevator deflection required for 1 g flight through the vor- 
tices is indicated. The limit and ultimate design load factors for the 
penetrating airplanes are indicated in figures 8, 9, and 10. 

The elevator deflection required for 1 g flight was defined as the 
deflection necessary to maintain the normal load factor at the 1 g datum 
value and is expressed by a relation developed from the second-order 
differential equation for the longitudinal motion of the penetrating air- 
plane. An example of the load factor induced by the vortices (elevator 
fixed) and the calculated elevator deflection required to maintain 1 g 
flight in this example are shown in figure 11. 
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Vertical Displacement 

The vertical displacements of airplane from prescribed penetra- 
tion paths were determined by numerically integrating, twice with respect 
to time, the normal acceleration (from 1 g flight) calculated for the 
airplane with elevator fixed. Vertical displacements of the penetrating 
airplane, as determined for several cases, are plotted in figure 12. 


DISCUSSION 


The discussion on the theoretical behavior of the trailing vortices 
is followed by accounts of the effects of certain variables on the normal 
load factors, the elevator requirement for 1 g flight through the vor- 
tices, and vertical displacements from the prescribed penetration paths. 


Trailing Vortices 

The vortices become essentially rolled up at specific distances 
behind the wings as listed in table IV. At these distances and beyond, 
the disturbances caused by the propeller and/or jet-engine-exhaust wakes 
were assumed to be negligible. The flight-test results of references 5 
and 6 indicate that this was a reasonable assumption. The magnitude and 
distribution of the resultant velocities of the trailing vortices as 
calculated in this report were found to be comparable to flight measure- 
ments given in reference 5* 

The flight-test results (refs. 5, 6, and 7) also show that the 
rolled-up vortices can maintain approximately initial strength for 30 sec- 
onds and up to possibly a minute or more in very calm air. No dissipation 
of vortex strength was considered because of the slow rate of decay of 
the vortices . Therefore the assumption was made that the spanwise wing 
loading was symmetrical for an airplane penetrating the vortices perpen- 
dicular to the vortex trail. Small amounts of atmospheric turbulence 
can, however, cause the vortices to move erratically and decay rapidly 
(refs. 5 and- 8). The circulatory velocity contours of the vortices were 
calculated from equations derived for wings having an elliptical spanwise 
loading. These equations would be applicable for wings having near 
elliptical loading either with or without flap deflections provided the 
flaps are full span. If the flaps are partial span, and deflected, the 
spanwise loading distributions would probably not be elliptical in shape. 
It is possible, then, that more than two rolled-up vortices could be 
developed rearward of each wing. 
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The behavior of the trailing vortices behind wings having abrupt 
changes in spanwise loading - that is, whether or not the trailing vor- 
tices will combine into two or more rolled-up vortices - may not be known 
(ref. 13). If knowledge of the trailing-vortex behavior is unknown, it 
is believed that a fair quantitative representation of the vortex-system 
velocity gradient can still be made, for use in determining the range of 
airplane normal load factors, by employing the equations presented herein. 


Load Factors 

Elevator-fixed penetrations .- The normal- load- fact or calculations, 
assuming rigid penetrating airplanes, include the load factors associated 
with the dynamic response of the airplanes to the vortex impulses . The 
assumption of a rigid structure for airplane appears reasonable 
(ref. 1^) for these calculations. However, if a penetrating airplane 
has a relatively flexible wing (such as, possibly, airplane Pg) the rigid- 
body treatment of the problem may not be adequate. Results of refer- 
ence 15 show that the wing-root bending moments and fuselage center-line 
load factors of an elastic-body system can be 25 percent greater than for 
a rigid-body system subjected to like flight conditions. The increases 
are attributed to a dynamic-overshoot effect as the airplane penetrates 
a velocity gradient . A rule of thumb is given in reference 15 for 
deciding whether or not the airplane, when encountering a sine-shaped 
gust, may be treated as a rigid body. The rule states that "as long as 
the period ratio, the time to penetrate the gust to the point of maximum 
gust velocity divided by one-fourth the natural period of the fundamental 
wing bending mode, is of the order of 5 or greater, the airplane may be 
treated as a rigid body." 

The normal load factors calculated for three penetrating-airplane 
velocities indicate two expected trends (fig. 7(a)). One trend is that 
peak load factors for a given airplane occur over corresponding shorter 
periods of time as the velocity of the penetrating airplane increases . 

The other trend is that the peak load factors show an increase in magni- 
tude with corresponding increases of airplane velocity. 

The load factors for the three center-of -gravity locations, illus- 
trating the effect of airplane static stability, show that the peak 
values of load factor vary little with center-of -gravity location 
(fig. 7 (b)). There are some differences in load factor between the 
peak values; however, the differences are of relatively small magnitude. 

The calculated load factors for different levels of penetration 
indicate that as the airplane becomes more distant from the line of 
vortex centers the peak load factors rapidly diminish in magnitude , as 
would be expected (fig. 7(c)) • Along penetration paths 25 feet above 


*\ 
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or below the line of vortex centers the peak load factors from the 1 g 
steady flight datum are less than one-half the peak values calculated 
for the path coinciding with the line of vortex centers . 

When the penetrating airplane enters the vortex system at angles 
to the line of vortex centers, and through the center of the core first 
approached, the peak maximum and minimum load factors experienced 
initially are approximately the same for entry angles to ±45°, the extent 
of the investigation. As the airplane progresses along the penetration 
path the load factors then are seen to approach the steady-state or 1 g 
value with increasing rapidity as the penetrating angle increases 
(fig. 7(d)). The load factors approach the steady-state value rapidly 
after the initial peak values since the influence of the farthest vortex 
is small for angles of ± 15 ° and ±45° . 

A comparison of the calculated load factors on the penetrating air- 
planes, traversing with elevator fixed the vortices of the three heavy 
transport airplanes (figs. 8 and 9)> indicates that the design limit and 
ultimate load factors can be, in some cases, exceeded. The design ulti- 
mate load factors are exceeded on airplane Pa when traversing the wakes 
of airplanes Gp and Gq and on airplane Pjj when traversing the wake of 
airplane Gc • The design limit load factor is exceeded on airplane Pg 
when traversing the wake of airplane Gg. The design load factors are 
not exceeded on either penetrating airplane when traversing the wake of 
airplane Ga- The peak load factors on airplane Pa are considerably 
greater than those on airplane Pp for comparable conditions because air- 
plane Pa is much lighter than airplane Pp. However, the ratios of the 
incremental peak load factors (from the 1 g datum) to the respective 
design ultimate load factors of the penetrating airplanes for comparable 
conditions are approximately the same. 

Estimated elevator motion .- The preceding discussion has been con- 
cerned with the penetrating airplanes traversing the vortex system with 
elevator fixed. Also investigated was the effect of an elevator control 
motion which might be made by a pilot in an effort to attenuate the 
normal load factors imposed by the vortices. As indicated in figure 10, 
the load factor n in this example exceeds both the design limit and 
the design ultimate load factors, whereas the design load factors were 
not exceeded when the same vortex system was traversed with elevator 
fixed. 

Of course an elevator motion could either aggravate or alleviate 
the normal load factor on the airplane traversing the vortices, depending 
upon the phasing of the elevator deflection with the load factors imposed 
by the vortices. However, since the trailing vortices are normally 
invisible and an encounter unexpected, a deflection of the elevators 
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would in a ll probability, because of inherent pilot and control-system 
lag, cause greater peak load factors than if the elevators were held 
fixed. 


Elevator Deflection Requirement 

It is also of interest to determine the elevator deflection necessary 
to maintain 1 g flight for the penetrating airplane while traversing the 
vortices. As an example (case 16) the incremental elevator deflection 
required to maintain 1 g flight was calculated for airplane P^. The cal- 
culations were made for case 1 6 primarily because the relatively small 
magnitude and low rate of change with time of airplane load factor for 
this case would require comparatively small incremental elevator deflec- 
tions and deflection rates . The low rate of change of load factor also 
permitted adequate solutions of n and ri (first and second numerical 
differentiations of n), which were needed to solve for the elevator 
deflections . 

As shown in figure 11, even if the required elevator motion could 
be anticipated by the pilot, the maximum deflection required to maintain 
1 g flight for the example case is not within the capabilities of the 
control system. 


Vertical Displacement 

The vertical displacements of airplane P^ from the penetration path 
prescribed for cases 1 to 5 range from approximately 4 feet above the 
line of vortex centers to 8 feet below (fig. 12). Specifically, the 
displacements for case 2, for example, range from about 3 feet above to 
8 feet below the line of vortex centers. From inspection of figure 7(c), 
which presents n for several values of h at the velocity of case 2, 
the effects of displacements of 3 feet or less in the region of the 
vortex cores would be expected to have little effect on the normal load 
factor. The 8-foot displacement which occurs well after the penetrating 
airplane passes through the second core would have negligible effect on 
the load factor. Thus it was concluded that, for the cases investigated, 
the effects of the vertical displacements on the normal load factors did 
not warrant inclusion in the present calculations . 


CONCLUSIONS 


The systematic calculation of normal load factor utilized available 
theories and accounted for normal translational velocity, lag in lift, 
and dynamic response in pitch of light airplanes traversing vortex fields. 
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The results indicate that light normal- cate gory and light transport- 
category airplanes traversing wakes of currently operational heavy trans- 
port airplanes can experience loading conditions that exceed the design 
limit and, in some cases, the design ultimate load factors. 

Load factors calculated for the light airplanes traversing the wake 
of a proposed supersonic transport were shown to exceed the design ulti- 
mate load factors by a comparatively large amount. Hence, because of 
the design trend toward heavy supersonic transports with short wing 
spans, penetrating airplanes can be expected to encounter loads of 
increased severity in the future . 

The results show that lighter airplanes can be expected to experience 
greater load factors than heavier airplanes when traversing a vortex sys- 
tem under similar conditions. The ratios of the maximum load factor (from 
level flight) to the design ultimate load factor were, however, about the 
same for the two penetrating airplanes used in this analysis. It was 
found, furthermore, that the load factors on the penetrating airplane 
increased with velocity, reached maximum values near the vortex cores, 
and were not affected appreciably by center-of- gravity location. It was 
shown that the penetrating airplane must pass either through or very close 
to the rather confined areas of vortex cores to experience load factors 
of severe magnitude. 

The results also show that although load factors imposed by the 
vortices could be alleviated by elevator deflection the pilot would be 
required to anticipate exactly the phasing of elevator motion with 
impulses from the vortices. However, because of inherent pilot and 
control-system lag, the loading conditions would almost certainly be 
aggravated rather than alleviated by elevator deflection. 


Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Field, Va., March 13 , 1961. 
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APPENDIX A 

SYMBOLS 

2 

a ^ normal acceleration, ft /sec 

A wing aspect ratio, b 2 /S 

b wing span, ft 

b-t horizontal-tail span, ft 

c wing mean aerodynamic chord, ft 

lift coefficient, nW/qS 

^L,t horizontal-tail lift coefficient, (L/qS)^ 

C m pitching-moment coefficient of airplane without horizontal 

tail, Mb/qS 2 

C m £ pitching-moment coefficient of isolated horizontal-tail 

surface 

d distance between core centers of the rolled-up trailing wing- 

tip vortices, ft 

Dp,D 2 ,D^,D 4 distances, ft (fig. 13 ) 

e distance behind wing at which trailing wing-tip vortices become 

essentially rolled up, ft 

F(t) impressed moment forcing function, M/-I, l/sec^ 

g acceleration due to gravity, 32.2 ft/sec^ 

h penetration-path height from line of vortex centers, ft 

I pitching moment of inertia of airplane, slug-ft^ 

ky radius of gyration about pitching axis, ft 

empirical constant denoting ratio of damping moment of complete 
airplane to damping moment of tail alone 


K' 



» . J 
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K]_,K 2 ^**.Kg constants occurring in equation (B2l) (defined in table V) 

L lift, lb 

m airplane mass, slugs 

M impressed airplane pitching moment about the center of gravity, 

ft-lb 

M w pitching moment attributed to circulatory flow of vortices, 

ft-lb (defined by eq. (B19)) 

Mg pitching moment attributed to elevator deflection, ft-lb 

(defined by eq. (B17)) 

n airplane normal load factor, a n /g 

n^ load factor for trimmed steady flight (n]_ = l) 

n r load factor due to airplane response 

n 1 load factor due to circulatory flow of vortices 

q dynamic pressure, lb/sq ft 

r radial distance from center of vortex core, ft 

r c radius of vortex core, ft 

S wing area, sq ft 

St horizontal-tail area, sq ft 

t time after airplane entry into vortex field, sec 

t]_,t 2 ,t^,ti(. times at which the airplane enters and departs the vortex 
cores, sec (defined by eqs. (B 8 ) and (B 9 )) 

V velocity of airplane, ft/sec 

V v normal translational velocity of airplane, ft/sec 

w resultant velocity of vortices perpendicular to penetration 

path, ft/sec 

w c tangential velocity of circulatory flow inside the vortex 

core, ft /sec 


wa downward velocity of rolled- up vortices , ft/sec 

w y tangential velocity of circulatory flow in vortex outside the 

core, ft/sec 

W airplane weight, lb 

x distance along the penetration path where w is approximately 

2 ft/sec, ft (fig. 13) 

x^. length from center of gravity of airplane to aerodynamic center 

of horizontal tail (positive when tail is rearward of center 
of gravity), ft 


a 

°t 

°e 

a l 

“r 

OOy 

Ow 

7 

AS 

e 

nt 

9 p, 02*65 


airplane angle of attack, radians 

horizontal-tail angle of attack, radians 

effective induced angle of attack, radians 

Otf modified for lag- in-lift effects, radians 

induced angle of attack due to airplane response, radians 

induced angle of attack due to nonnal translational velocity 
of airplane , radians 

induced angle of attack due to resultant velocity w of 
vortices, radians 

angle of penetration path with respect to line of vortex 
centers, deg 

incremental elevator angle, radians 

j _ 

downwash angle, — ■ - a, radians 
da 

tail efficiency factor, /q. 
angles, deg (fig. 13) 

effective airplane pitch angle due to vortex circulatory flow, 
radians 


p mass density of air, slugs/cu ft 

The notations a, 0, h, n, and a denote first and second deriva- 
tives of variables with respect to time. 



Subscripts : 

g generating airplane 

p penetrating airplane 


t 


horizontal tail 


I 


16 


APPENDIX B 

DETAILS OF METHODS USED 


To formulate the methods used in this paper, well-known theoretical 
descriptions of wing-tip trailing vortices, as developed by a generating 
airplane, are first resolved to convenient form. The angles of attack on 
the penetrating airplane induced by the vortex flow were then determined 
and used to calculate the normal load factors . 


Trailing Vortices 


The transport airplane generating the vortices was assumed to have 
elliptical spanwise wing loading. In accordance with the following 
expressions from theory (ref. 12), each wing-tip trailing vortex was 
essentially rolled up at a distance e behind the wing, where 


e 


Ag 

0.280 b c 

Cl 8 


(Bl) 


and had at the initial rolled-up condition a strength (tangential velocity 
gradient of the circulatory flow) w v such that 


w v = 



(r ^ r c ) (B2) 


The tangential velocity of circulatory flow within the core of each 
vortex (r = r c ) for the initial rolled-up condition was considered to vary 
linearly from a maximum at the periphery (r = r Q ) to zero at the core 
center (r =0) as 


_ ( 2L \ r 

Wc ( 2 ) r c 

Vt bpr c V/g c 


The core radius is defined by the theory as 


(B3) 


' °' 155 ( 2 ) 


(B4) 



T 
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and the distance between core centers of the vortices is given by 




(B5) 


The vortex cores have a downward velocity with respect to the flight 
path as calculated from the relation 


w d 


-4 /Cp\ 

” Tt 2 V§ \ rtA /g 


(b6 ) 


A method was set up to evaluate the resultant velocities, contributed 
by the circulatory flow of the two vortex fields, perpendicular to the 
penetration paths. The resultant velocities at given times along the 
penetration paths parallel to the line of vortex centers outside the vor- 
tex cores (h ^ r c ) can be expressed by 




2L 


rt 2 bpV j 


[cos 0- 


2 \ 

cos Gg’ 

d 2 y 


(B7) 


where 


D 1 = x - Vpt 
D 2 = D 1 + d 

. , _i h 

0]_ = tan — 


0 2 = tan"^ -~ 


as shown in figure 13(a). Also for the cases where the penetration paths 
coincide with the line of vortex centers (h = 0) the resultant velocities 
were calculated with equation (B7) except those within the cores. The 
resultant velocities within the cores were defined by the straight line 
joining the resultant velocities, calculated with equation (B7), at the 
times the penetrating airplane entered and departed the cores. The times 
the first core was entered and departed are given, respectively, as 


* _ x - r c 

1 = _ V 

x + r P 


and 


(B8) 
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and for the second core the times are 


and 


x + d - r c 



v 4- rl 4- r» 



(B9) 


The resultant velocities perpendicular to penetration paths at 
angles to the line of vortex centers about the vortex center first 
approached by the penetrating airplane are, except within the vortex 
core , 


where : 





x - v p t 


D, = Di + d cos 7 


= d sin 7 
9 




= tan 1 _ 


(BIO) 


as shown in figure 13(b). Resultant velocities in the region of the core 
entered were also defined by the straight line joining the resultant 
velocities, calculated with equation (BIO), at the times the airplane 
entered and departed the vortex core. The equations used for expressing 
the times to enter and depart the core, tq and t 2 > respectively, are 
the same as previously written (eqs. (B 8 )). In this analysis the penetra- 
tion paths at angles to the reference line do not intercept the second 
core; hence t^ and t 4 are not considered. 

An example of the velocities w calculated with equations (B7), 

(B 8 ), and (B9) for a penetration path coinciding with the line of vortex 
centers is given in figure 5 (case 2). In the subsequent calculations the 
velocities w are assumed to be of uniform strength across the wing span 
of the penetrating airplane at any instant of time during penetration of 
the vortices. 



1 
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Equations (B7) and (BIO), then, in conjunction with equations (B8) 
and (B9) , are used in this report to describe the velocity gradients 
experienced by the penetrating airplane while traversing the trailing 
wing-tip vortices. The effect of the vortex-flow resultant velocity com- 
ponent parallel to the penetration paths on the load factors of the pene- 
trating airplane was believed to be small and was therefore not included 
in the calculations. Also, the interaction between the vortices of the 
generating and penetrating airplanes was considered negligible. 


Normal Load Factors 

The normal load factor n on the penetrating airplane, while tra- 
versing the vortex fields, can be written as 


n = n^_ + n^ + n r 


(Bll) 


and in coefficient form as 


n ■ 1 + (sr t) p (a * + (B12) 

To solve equation (B12) for n the incremental induced angles (a i and 
Oy) must first be dete mined. The values of the other parameters in the 
equation are known for the penetrating airplane for which calculations 
are to be made. 

Induced angle of attack ctq.- Two steps were taken to resolve the 
resultant vortex velocities w to the induced angle of attack ctq on 
the penetrating airplane (no pitch motion assumed) while traversing the 
trailing vortices . The first step taken evaluated 



at times 


t = 


x - D]_ 


V P 


(B13) 


(B14) 


along the penetration path. Several examples of are shown plotted 

in figure 14 for the different types of cases investigated. The 



r 
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penetrating airplane was considered to enter the vortices at time t = 0 
(arbitrarily selected at the point where w ^ 2 ft/sec). 

The second step evaluated a j by making adjustments for the lag- 
in-lift effects (Kiissner function) on the induced angle ay. The pro- 
cedure utilized was a superposition integral that approximated Oy by 
a series of step changes in angle of attack. The function of figure 15 
(obtained from ref. 16) was applied to each unit change to account for 
the lag in lift. As illustrated in the figure, the angle -of -attack ratio 
a z/ a w becomes unity at 8 wing-chord lengths of penetration into a step- 
velocity gradient. Examples of the induced angles of attack and a 2 

on the penetrating airplane are shown in figure 16 for case 2. Also shown 
in the figure is the effective induced angle of attack a e which will be 
derived later. 

Induced angle of attack Oj. . - The induced angle of attack cxj. can 

be expressed in the differential equations for the longitudinal motion of 
an airplane, namely, equations summing the forces perpendicular to the 
relative wind and the pitching moments about the airplane center of 
gravity. Implicit in these equations are the assumptions of a rigid air- 
plane, linear aerodynamic derivatives, and small changes in penetrating- 
airplane speed and penetration-path angles while traversing the vortices 
along a specified path. The equations are reducible (ref. 17) to the 
equivalent second-order differential equation of longitudinal motion 

a*. + Kqoj. + K 20 r + = 0 (B15) 

The term M is the impressed pitching moment on the airplane about 
the center of gravity; Kq and Kp are constants for any given set of 
flight conditions as defined in table V. 

The impressed pitching moment on the penetrating airplane under the 
present considerations is the sum of the moments attributed to the circu- 
latory flow of the vortices My. and the moments due to elevator deflec- 
tion Mg. Hence 


;j = (>% + m»> h 


(B16) 
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where 


— - -AS £0- d ° L > t „ SiS. + dC m,t S t 2 

-I ' 2m ^ d6 ''t dB V bt% 2 

dc L,t dd L,t K'lt 2 p 

da-j (15 ^2 mky 2 / 

which can be written in abbreviated form as 


and 


-i 


Mw s 

= oe - py2 

d c m : 

3 2 

-I 

^ 2m 

^- a bky^ 


• PV 

dC L,t 



0,(5 2m 

da t 



+ 9 e — 
e 2m 

dC L,t 

da t 

\ 


-K-jA6 


dC L,t S t x t / n de\ 

da t Tlt k Y 2 \ da / 

^t x t 2 de 
k Y 2 da 

S t x t^ K’ 

k Y 2 f ' it 


which in abbreviated form is written 


~ = K 4 a e + K 5 d e + K 6 9 e 


(BIT) 


(BIB) 


(B19) 


(B20) 


The positive directions of the sign conventions employed are shown in 
figure 17 and are the same as those used in reference 17. The K con- 
stants , as defined, are listed in table V. 

Equation (B17) represents the moment about the airplane center of 
gravity attributed to an incremental elevator deflection from the trim 
position. In equation (B19) the term containing Og represents the 
moment on the airplane caused by the effective angle of attack induced 
by the resultant velocity gradient of the vortices; the term was 

introduced to represent the moment due to lag in downwash at the tail; 
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and the 6 e term represents the moment due to change in horizontal-tail 
angle of attack due to effective airplane pitch. 

Equation (B15) can he rewritten, using equations (Bl6), (Bl8), 
and (B20) as 

cLj. + Kqd-j. + = KjAS - (Ki^a^ + + K^0 e ) (B2l) 


which is the longitudinal equation for a damped pitching oscillation 
with the impressed moment forcing function 


F(t) = K 5 AS - (K4ot e + K 5 d e + Kg0 e ) 


(B22) 


The longitudinal equation of motion (eq. (B2l)) is solved for the 
induced angle of attack Op, the response of the airplane about the 
penetration path while traversing the vortex system. To solve the equa- 
tion, the evaluation of the K constants and time histories of the 
effective induced angle of attack a e (no pitch motion), de> 6 e , and 
incremental elevator deflection A£>, if any, are required. The values 
of the K constants as determined for the flight conditions investigated 
are listed in table VI. 

To resolve the previously determined induced angle ctq to the 
effective induced angle of attack a e , an adjustment was made for the 
effect of the penetrating airplane's normal translational velocity on 
the induced angle of attack aq (no pitch motion). The airplane's 
normal acceleration along the penetration path, attributed to the vortex 
velocity gradient, was calculated from 


% = ng 


'<1Cl qS \ 
do, W / 


oqg 


(B23) 


The normal translational velocity V y was then derived by integrating 
the acceleration a with respect to time, and the relation 


Ov = 



(B2*0 


was used to express the induced angle of attack on the penetrating air- 
plane due to its normal translational velocity. 
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Hence, the effective induced angle of attack on the airplane while 
passing through the vortices with no pitch motion can be written 


a e = a Z + o-v 


(B25) 


Time histories of the effective induced angle of attack a e were deter- 
mined for each case investigated. (See fig. 16 for a typical example.) 

The methods used to evaluate the rates of change of effective angle 
of attack dg and effective pitching velocity 8 g induced on the pene- 
trating airplane by the angle of attack a„ are illustrated in figure 18. 
The values of ag and 9 e were calculated from the expressions 

d e = (B26) 

x t 

and 

4 



Sufficient information is now available concerning the impressed 
moment forcing function F(t) (eq. (B22)) to solve the longitudinal 
equation of motion (eq. (B2l)) for the response Op of the penetrating 
airplane while traversing the vortex system. 

An example of a typical elevator-fixed impressed moment forcing 
function F(t) calculated with equation (B22) for case 2 is shown in 
figure 19. The forcing-function time histories of cases 1 to 16 were 
similarly found, and the time histories were programed according to 
equation (B2l) in an automatic digital computing machine. Initial con- 
ditions specified for equation (B2l) were dp = 0, dp = 0, and dp = 0 
at t = 0. Solutions to equation (B2l) for ctp were obtained at times 
selected along the penetration paths to represent adequately the airplane 
response. The unsteady-lift effects following a change in angle of attack 
(Wagner function) were relatively small and not included. 

General solutions to equation (B2l) are presented in chart form 
(fig. 20) for a range of airplane characteristics, and an example 
illustrating their use is given in appendix C . 


Elevator Deflection Requirement 


The incremental elevator deflection required for the attenuation of 
the load factor to 1 g flight through the vortex system is expressed by 
the relation given in reference 17 (with negative sign affixed) as 


AS = 



+ K x n + K 2 (n - l) 


(B28) 


Sufficient information about the terms in equation (B28) is already 
available to solve for the incremental elevator deflection required of 
the penetrating airplane. The load factor n on the penetrating air- 
plane traversing the vortices with elevator fixed is determined from 
equation (Bll). The n and h terms are determined by taking first 
and second derivatives of n. 


Vertical Displacement 

The vertical displacement of the airplane from the prescribed pene- 
tration path was obtained as a first approximation by integrating (numer' 
ically) twice with respect to time the normal acceleration a n on the 
penetrating airplane where a n = (n - l)g. 
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APPENDIX C 

EXAMPLE OF CHART APPLICATION 


An example is given in this appendix to illustrate the generality 
of the charts in figure 20 for computing the normal load factor on the 
penetrating airplane P A responding to an elevator deflection. The charts 
are solutions of equation (B2l) , in terms of airplane angle-of-attack 
response with unit forcing function l(t) (that is, F(t) = 0 

at t < 0 and F(t) =1 at t = 0). The ranges of and K 2 repre- 

sent values of airplane aerodynamic and geometric characteristics and 
flight conditions that are likely to exist. 

The normal load factor n r is defined as 

nr = (lr §5 ) p K20t (C1) 

where K 2 a r is defined in the relation 


’dj. + K^ctj. + K^Oj. = KjAS - (Kkc^ + + Kg0 e ) (C2) 


Physical and aerodynamic characteristics of the airplane used in the 
example are given in table II and the flight conditions are the same as 
those of case 2. (See table III.) 

The values of K^_ and as defined in table V, are also the same 

as those of case 2 . (See table VI.) Interpolating in figure 20(e) 

(Kx « 6.0 and K 2 « 21 . 6 ) gives the response KpOj. of the airplane to 
a unit step function. This response is shown plotted in figure 2l(a). 

The next step is to evaluate the impressed moment forcing function 
F(t) on the airplane from 


F(t) = KjAS - (K 40 ^ + K^ctg + Kg9 e ) 


(C3) 


The impressed moment used in this example is due only to an estimated 
incremental elevator deflection (shown plotted in fig. 10 (a)); hence, 
a-e = 0, 0 ^= 0 , and 0 e = 0. From table VI, = -32.5214-3. The 
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time history of the derived impressed moment forcing function is plotted 
in figure 2l(b). 

Unit step functions are then superimposed on the plot of derived 
impressed moment forcing function so as to represent adequately the 
curve. At selected times the responses KpOr of each unit step are 
algebraically summed and the sums are each multiplied by 0 . 5111 , the 

value of /— i for the example, to obtain the load factor on the 

\da K 2 w; p 

airplane n r due to the impressed moment forcing function. 

A comparison is made in figure 2l(c) between the airplane normal 
load factors n r calculated with the automatic computing machine and 
those determined by the superposition method using the charts. The 
agreement is considered to be good and illustrates that the charts are 
adequate for determining the response and load factors of an airplane 
with impressed moment. In a similar manner the incremental normal load 
factors could be calculated for other airplanes and/or flight conditions 
with the use of these charts. 
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TABLE II 

PEI^TRATIKG-AJRPLANE CHARACTERISTICS 

Airplane Airplane 

Geometric ?A 


Airplane weight, W, lb 2,000 35*000 

Gross wing area, S, sq ft 147*5 754.0 

Wing span, b, ft 29-3 95-2 

Wing aspect ratio, A 5-82 12.00 

Wing mean aerodynamic chord, c, ft 5*25 8.45 

Horizontal tail area, S^, sq ft 27.2 172.0 

Horizontal tail span, b-j., ft 9*5^- 31*65 

Tail length, x^, ft: 

Center of gravity, 20 percent c 12.11 

Center of gravity, 25 percent c 11.85 

Center of gravity, 30 percent c 11*59 35*70 

Radius of gyration, ky, ft 4.22 15*50 

Elevator deflection, 5, radians: 

Maximum up 0.4188 

Maximum down 0 . 2269 

Aerodynamic 

Slope of airplane lift curve, dC^/da, 

per radian 4.18 5*17 

Slope of tail lift curve, dCj, f/dou-, 

per radian . 3*66 4.32 

Downwash factor, de/da 0.445 0.352 

Tail efficiency factor, 1.00 1.00 

Empirical airplane damping factor, K' 1.10 1.10 

Elevator effectiveness factor, dCp, ^/db, 

per radian 9 2.93 1*73 

Rate of change of tail moment with camber due to 

elevator angle, dC m ^-fc/db* per radian ....... -0.57 - 0.50 

Rate of change of moment coefficient with angle 
of attack for airplane without tail, dCj^/da, 
per radian: 

Center of gravity, 20 percent c 0.075 

Center of gravity, 25 percent c 0.284 

Center of gravity, 3° percent c 0.490 1.118 


PENETRATING-AIRPLANE FLIGHT CONDITIONS 
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TABLE VI 

VALUES OF K CONSTANTS OCCURRING IN LONGITUDINAL EQUATIONS 
OF MOTION FOR THE VARIOUS CASES 


Case 

*1 

k 2 

k 3 

K4 

k 5 

k 6 

1 

3.96273 

9.52813 

-14.37206 

6.30771 

0.75701 

1.87150 

2 

5.96240 

21.56047 

-32.52143 

14.27322 

1.13902 

2.81590 

3 

6.94393 

29.25070 

-44.12123 

19-36422 

1.32652 

3.27945 

4 

5.79076 

15.85335 

-31.82337 

8.88237 

1.08958 

2.69369 

5 

4.60787 

27.34639 

-33-22020 

19.73585 

1.18955 

2.94082 

6 

5.96240 

21.56047 

-32.52143 

14.27322 

1.13902 

2.81590 

7 

5.96240 

21.56047 

-32.52143 

14.27322 

1.13902 

2.81590 

8 

6.9^393 

29.25070 

-44.12123 

19.36422 

I.32652 

3.27945 

9 

6.94393 

29.25070 

-44.12123 

19.36422 

1.32652 

3.27945 

10 

6.94393 

29.25070 

-44.12123 

19.36422 

1.32652 

3.27945 

11 

6.94393 

29.25070 

-44.12123 

19.36422 

1.32652 

3-27945 

12 

3.45413 

7.48648 

-5.76739 

5.30838 

.49809 

1.55654 

13 

3.45413 

7.48648 

-5.76739 

5.30838 

.49809 

1.55654 

14 

3.45413 

7.48648 

-5.76739 

5.30838 

.49809 

1.55654 

15 

5.96240 

21.56047 

-32.52143 

14.27322 

1.13902 

2.81590 

16 

3.96273 

9- 52813 

-14.37206 

6.30771 

•75701 

1,87150 
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Figure k.- Example of initial circulatory velocity w v contours of rolled-up trailing vortices 
in the plane of penetration paths. Calculated for airplane G A . 
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(a) Variation of penetration velocity. 



Time, sec 


(b) Variation of center-of -gravity location. 



Time, sec 


(c) Variation of penetration level. 



Time, sec 

(d) Variation of penetration angle. 

Figure 7*- The effects of several variables on the calculated load 
factors on airplane with elevator fixed traversing the rolled 
up trailing vortices of airplane Ga. 







Load factor 



Time, sec 


Figure 8.- Calculated, load, factors on airplane with elevator fixed 
traversing the trailing vortices of airplanes Ga, Gg, and G^ . 
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Figure 9-- Calculated load factors on airplane Pp with elevator fixed 
traversing the trailing vortices of airplanes Ga, Gg, and Gq. 
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(a) Penetration paths parallel to the line of 
vortex centers at a distance h. 



(b) Penetration paths at an angle 7 with 
the line of vortex centers. 

Figure 13- - Conventions used in calculating the resultant vortex veloci- 
ties w perpendicular to the penetration paths. Positive directions 
of quantities shown. 
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(a) Flight path on line of vortex centers. 



(b) Flight path at a distance parallel to the line of 
vortex centers, h = ±15 ft. 



(c) Flight path at an angle to the line of 
vortex centers. 7 = ±15° • 


igure Examples of angle of attack induced on penetrating airplane 

by the trailing vortices for the different types of cases investigated, 


176 ft/ 






Figure 15 •- Development of induced angle of attack on a wing penetrating a sharp-edge gust (lag 

in lift). A = 6. (From ref. 1 6 .) 




Figure l 6 .- An example of the effects of lag in lift and airplane normal translational velocity 
on the induced angle of attack on the penetrating airplane. 0^=0^ -a 1 ; case 2 . 
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Figure 19.- Example of a typical elevator-fixed impressed moment forcing function F(t). Calcu- 
lated for case 2: AS =0. 












Figure 20.- Concluded. 





(b) Superposition of F(t) with unit step functions l(t). 



(c) Comparison of calculated airplane load factors. 

Figure 21.- Example of applying the derived charts of figure 20 to calcu- 
late a load factor due to F(t). 
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